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A simple empirical model and approach are introduced for calculation of the vibrational spectra 
of arbitrary single wall carbon nanotubes. Differently from the frequently used force constants 
description, the model employs only invariant quantities such as variations of lengths and angles. 
All the salient qualitative features of vibrational spectra of nanotubes naturally follow from the 
vibrational Hamiltonian of graphene upon its isometric mapping onto a cylindrical surface and 
without any ad hoc corrections. A qualitative difference with previous results is found in a parabolic, 
rather than a linear, long wavelength dispersion of the transverse acoustic modes of the nanotubes. 
The parabolic dispersion is confirmed and elucidated in the provided continuum analysis of the 
vibrations. We also discuss and use an alternative definition of the nanotube unit cell with only 
two carbons per cell that illustrates a "true" longitudinal periodicity of the nanotubes, and of the 
corresponding Brillouin zone. 

PACS numbers; 61.46.+W, 62.25.+g, 62.30.+d, 46.40.-f 
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I. INTRODUCTION 

Vibrational spectra of individual single wall carbon 
nanotubes are of considerable interest and have been cal- 
culated previously within different frameworks such as 
an empirical force constant model, ah initio studies'* 
and tight-binding molecular dynamics.^ It is known that 
the higher-frequency part of the nanotube spectra is rel- 
atively well represented already by the zone folding of 
the graphene spectrum. The lower-frequency part, how- 
ever, has generic features owing to the one-dimensional 
character of nanotubes. Particularly, the spectra exhibit 
four types of acoustic modes with vanishing frequencies: 
one longitudinal, two transverse and one twisting. The 
existence of these modes has to do with general consider- 
ations ~ displacements of the tube as a whole along and 
perpendicular to its axis, and the rotation of the tube 
about the axis do not cost energy - rather than with spe- 
cific nanotube interactions. Analogous vibrations were 
also discussed in the context of quantum wires as dilata- 
tional, flexural and torsional modes. ^ These modes are 
important contributors to the low-temperature quantized 
thermal conductance of such phonon waveguides. ^'^ 

The frequently used force constant model of Ref. 1 was 
first developed for planar graphene based on the experi- 
mental data for graphite and then adapted for nanotube 
geometries. Direct application of the graphene force con- 
stant values was not found to lead to zero frequencies 
for all four modes mentioned above. To overcome this 
difficulty, special curvature corrections were introduced 
to the force constants.* In this paper we develop another 
empirical model, where the harmonic vibrational Hamil- 
tonian of graphene is built using only "invariant" quan- 
tities such as variations of bond lengths, interbond and 
dihedral angles. Such a description is similar in spirit to 
used in conformational analysis and stereochemistry (see, 
e.g., Ref. 9) and in bond models for vibrations in covalent 
semiconductors (Ref. 10 and references therein). Further 
mapping of the graphene Hamiltonian onto a cylindrical 



surface of the nanotubes then allows to derive all fea- 
tures of arbitrary nanotube spectra "naturally" , without 
any curvature corrections. This way the idea of Ref. 1 of 
using the same type of vibrational Hamiltonian for both 
graphene and nanotubes turns out to be realized with no 
need for ad hoc modifications. 

We find a qualitative difference with the previously 
published results in a parabolic, rather than a linear, dis- 
persion of the transverse acoustic modes of nanotubes. 
This parabolic dispersion is further illustrated in the 
corresponding continuum model of vibrations and is in 
agreement with the analysis of vibrations of elastic cylin- 
drical shells.** The continuum model also shows the ori- 
gin of another salient long wavelength feature of the nan- 
otube spectra: a coupling of the longitudinal acoustic 
with the breathing mode. We believe the parabolic char- 
acter of the low-frequency part of the transverse mode 
dispersion is quite generic similarly to the well-known 
bending waves of rods*^ with wavelengths much longer 
than the rod size. Such a parabolic dispersion was, 
e.g., calculated for the lowest flexural modes of quan- 
tum wires;^'^ recent applications of elastic cylinder mod- 
els also include vibrations of cytoskeletal filaments and 
microtubules.*^"*^ To our knowledge, the parabolic dis- 
persion was not calculated previously for carbon nan- 
otubes. 

With our approach we can easily calculate vibrational 
spectra of arbitrary (TV, M) nanotubes. In doing this, we 
also employ an alternative definition of the nanotube unit 
cell with only two carbons per cell as opposed to possibly 
many carbon atoms of the conventional definition.* This 
way a "true" longitudinal periodicity of the nanotubes is 
elucidated. The period is a projection of one of the prim- 
itive vectors onto the nanotube axis, the primitive vector 
itself being in general not parallel to the axis. The Bril- 
louin zones can correspondingly be wider than usually 
used and the total number of the vibrational branches 
turns out to be just 6dc, where dc is the greatest com- 
mon divisor of N and M . 
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For numerical computations, we will be using some 
parameterization of the model, specifically based on a 
set of data from Ref. 1 as well as on some experimental 
data. However, the numerical computations here serve 
mostly illustrative and qualitative purposes. 

II. VIBRATIONAL HAMILTONIAN 

Despite the fact that atoms of carbon nanotubes are 
arranged in a 3-d fashion, the excitations of nanotubes 
can be described in the way very similar to excitations 
of the planar graphene. To clearly see this connection 
for vibrational excitations, one can use local (position- 
dependent) coordinate systems for atomic displacements 
and the clastic potential energy written in an "invariant" 
form. This way the qualitative transformation from vi- 
brational spectra of an infinite plane to spectra of curved 
cylindrical structures appear naturally without ad hoc 
corrections. 

Following Ref. 1, we also consider carbon-carbon elas- 
tic interactions up to the fourth nearest neighbor. The 
harmonic potential energy U = Ui + Uo is however ex- 
pressed as a function of only invariant quantities such as 
variations of bond lengths and various angles. The first 
term Ui would correspond here to in-plane deformations 
of graphene and in general requires ten elastic parameters 

(ijk) 

+ K^SifiijkiSkj + Sljk)] (la) 

4 {il}—m 
m=3 (ijkl) 

+ K^iS(pijkSlki + SipjkiShj)] , (lb) 

while term Uo would describe out-of-plane graphene dis- 
tortions and needs three elastic parameters: 

4 {il}—m 

^°=E E ^TSxlki- (2) 

m=2 (ijkl) 

The structure of Eqs. (1,2) can be conveniently thought 
of in terms of triangular plaquettes (ijk) formed by bonds 
(ij) and (jk) connecting nearest carbons i and j, and 
j and k, respectively. Variation of bond (ij) length 
is denoted Slij , and variation of the inter-bond angle 
at the common carbon j denoted Sipijk- Correspond- 
ingly, Eq. (la) completely describes the deformation en- 
ergy of individual plaquettes. Equations (lb, 2), on the 
other hand, completely describe the interactions of de- 
formations on neighboring plaquettes; specifically, (ijkl) 
stands for plaquettes (ijk) and (jkl) that are adjacent 
along bond (jk). There are three different ways to form 
adjacent plaquette pairs and notation {il} = m distin- 
guishes them by indicating that carbons i and I are the 




FIG. 1: Nomenclature of possible triangular plaquette pairs 
adjacent along a nearest-neighbor carbon-carbon bond. To 
guide the eye, individual adjacent plaquettes are filled with 
lines at different angles. 



mth nearest neighbors (m ~ 2,3,4), as illustrated in 
Fig. 1. (In case to = 2, it is actually plaquettes (ijk) 
and (Ijk) that are adjacent along (jk). This case is not 
explicitly included in Eq. (lb) because of the constraint 
that a sum of inter-bond angles for three plaquettes sur- 
rounding a carbon atom is fixed.) The "out-of-plane" 
interaction, Eq. (2), involves the variations Sxijki of the 
dihedral angles between corresponding plaquettes. 

It is worth stressing that for the planar graphene, 
Eqs. (1,2) give the most general description of harmonic 
interactions involving up to the forth nearest neighbors. 
As such, it, of course, can reproduce a force constants 
description. The latter would be derived by expanding 
Eqs. (1,2) in the atomic pair differences. Evidently, then 
there would be certain relationships between a larger 
number of force constants as dictated by the invariance of 
the potential energy with respect to overall rotations.^® 
Specifically, ten elastic parameters of Eq. (1) yield twelve 
force constants, comprising eight in-plane constants of 
the type explicitly considered in Ref. 1 and four con- 
stants mixing radial and tangential displacements that 
were implicitly set to zero in that reference. Three elas- 
tic parameters of Eq. (2) yield four out-of-plane force 
constants.^ 

The advantage of using invariant quantities in 
Eqs. (1,2) is that the same functional form of the po- 
tential energy can be directly used when carbon atom 
positions are (isometrically) mapped from the graphene 
plane onto the cylindrical surface of a nanotube. When 
on the nanotube surface, variations of the bond lengths 
and various angles just need to be calculated from the 
carbon atom displacements using the actual curved ge- 
ometry. Of course, a new set of force constants appro- 
priate for the now curved geometry can be again derived 
through the pair expansions. Force constants so obtained 
would automatically obey the correct relationships to sat- 
isfy the invariance with respect to rotations. 

The standard translational invariance of the vibra- 
tional Hamiltonian is preserved if displacements of car- 
bons are expressed not in terms of common (xyz) co- 
ordinates but in terms of local orthogonal coordinates 
(uvw): u - along the nanotube axis, v - perpendicular 
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to the tube axis and parallel to the tube surface, and w 
- perpendicular to the tube surface. These simple ideas 
are illustrated in more detail in Appendix A for an easier 
to follow example of the relationship between vibrations 
of a linear chain of atoms and vibrations of a ring of 
atoms. For the problem at hand, we use the local dis- 
placements {uncVnajWna) depending on carbon na ac- 
tual geometric position on the cylindrical surface of the 
tube. (Here carbon index na consists of a 2-d vector n 
specifying the unit cell of the parent graphene plane and 
a = 1,2 specifying one of the 2 carbons in the graphene 
unit cell.) With local displacement bases in place, the 
invariance of the vibrational Hamiltonian with respect 
to translations by graphene primitive vectors is held on 
equal footing in graphene and nanotubes: wave vectors k 
will "know" only differences between neighboring carbon 
indices na. The problem is thereby reduced to calcula- 
tion of the usual, "graphene-like" ,6x6 dynamical ma- 
trix, but which naturally contains the correct mixing of 
the "in-plane" and "out-of-plane" displacements in nan- 
otubes. Using proper quantization rules for the allowed 
phonon wavevectors, one can then readily derive the vi- 
brational spectra of nanotubes of arbitrary chirality. 

III. PRIMITIVE CELLS AND BRILLOUIN 
ZONES 

The geometry of a single wall carbon nanotube is de- 
termined by the chiral vector^ 

Ch = Nai+MaL2, (3) 

where ai and a2 (|ai| = |a2| = a) are two prim- 
itive vectors of the 2-d graphene (hexagonal) crystal 
structure. Vector C/i is perpendicular to the nanotube 
axis. Conventionally,^ the translational vector parallel 
to the tube axis is defined, which we denote here as 
Td: Td = hai + t2ai, where ti = (2M + N)/dR, 
t2 = — {2N + M) jdji and rf/j is the greatest common di- 
visor of (27V -I- M) and (2M-I-7V). The resulting unit cell 
of the nanotube built of C/j and T/5 can contain many 
carbons Na = 4:{N^ + M"^ + NM)/dR, and the longitu- 
dinal period \Td\ of chiral tubes be much larger than a. 
The corresponding Brillouin zones (BZs) would then be 
narrow and contain many excitation spectrum (whether 
vibrational or electronic) branches. 

In this paper we use an alternative picture of the unit 
cell and BZ construction that is aimed at having as small 
number of branches in the zone as possible. As described 
in more detail in Appendix B, this number of branches is 
determined by the greatest common divisor of N and M, 
denoted by dc- The total number of continuous branches 
in BZ is equal to 2dc per each degree of freedom of a car- 
bon atom, that is, 6dc for vibrational excitations. This 
corresponds to the nanotube unit cell containing only two 
carbons and which can, e.g., be built with primitive vec- 
tors Cfi/dc and T. Different from vector Td, the trans- 
lational vector T = Pai -|- Qa.2 is in general not paral- 



lel to the nanotube axis. As discussed in Appendix B 
(see Eq. (B7)), integers P and Q here satisfy condition 
MP — NQ = dc- It is the projection of vector T onto 
the axis that determines the longitudinal period 

0,1 = sfld^c?l2Ch (4) 

and the width 27r/a[| of the BZ. 

The number of branches in the BZ is related to the 
transverse quantization of the 2-d wave vector k of the 
parent graphene band excitations: 

kCft = 21,1, (5) 

resulting in the appearance of one-dimensional sub-bands 
characterized by the integer quantum number I. The con- 
struction of the unit cell and BZ employed in this paper 
recognizes that there would be only d^ "unique" quan- 
tization levels (that is, only dc independent integers /): 
all other allowed k-vectors can be obtained with trans- 
lations by graphene reciprocal vectors. Many branches 
of the conventional^ BZ would not exhibit gaps at that 
BZ boundary; they would correspondingly become sin- 
gle continuous bands when properly "unfolded" in our 
construction. The quantization lines within our BZs can 
span several hexagons of the graphene reciprocal lattice 
(see example of Fig. 9). 

IV. VIBRATIONAL SPECTRA 

As was mentioned above, potential energy in Eqs. (1,2) 
is capable of reproducing results of the force constant 
model^ for the graphene spectrum. Figure 2 (a), how- 
ever, has been calculated with a parameterization of elas- 
tic constants Ki in Eqs. (1,2) such as to achieve only a 
close similarity to the published spectrum.^ In our cal- 
culations, we chose to slightly and somewhat arbitrar- 

ily modify the tangential force constants from the 
published^ values so as to satisfy -\- 6(^j^^ -|- 4^^' -|- 
14(/)f^^ = 0. The latter equality is required by the ro- 
tational invariance - elastic energy should be zero for 
the overall rotation of the graphene plane. Original 
constants^ do not obey it. A recently published new set 
of force constants^'' also does not satisfy this require- 
ment. The overall scaling of elastic constants was chosen 
here so as to reproduce graphene experimental optical 
frequencies of 1580 and 868 cm^^. Once the values of 
elastic constants Ki have been defined for graphene, the 
same values are used to calculate spectra of nanotubes, 
examples of which are shown in Figures 2 (b)-(d). As 
discussed in Sec. II, we do not need the knowledge of 
force constants of Ref. 1 for curved geometries because 
the Hamiltonian used automatically preserves all invari- 
ance requirements. 

In displaying the nanotube spectra, we use the defi- 
nition of BZ as discussed in Sec. HI and Appendix B. 
With that definition, the BZ contains only 6(ic vibra- 
tional branches. For armchair tubes, ay = a/2 and the 
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FIG. 2: Model vibrational spectra: (a) Graphene; (b) (10,10) 
tube; (c) (10,9) tube; (d) (16,1) tube. See text for definition 
of ail. 



BZ in Figure 2 (b) is twice as wide as the conventional 
BZ, that picture would be restored by simple folding. 
Figure 2 (b) exhibits dc — 10 "independent" transverse 
quantization levels. On the other hand, Figures 2 (c) and 
(d) with dc = 1 correspond to only one transverse quan- 
tization level, that is, to six branches. Each polarization 
band/branch exhibits a continuous evolution of the po- 
larization vectors with in general a strong /c|| dependence. 
It is worth noting that transverse acoustic modes in this 
picture have their frequency vanishing at finite (which, 
of course, lie in a hexagon of the graphene reciprocal lat- 
tice other than the one where fcy = is). 

Figure 3 shows low-frequency parts of the vibrational 
spectra of the (10,10) and (10,0) tubes reduced to a sin- 
gle hexagon of the reciprocal lattice. Panel (a) can be 
directly compared to the published results^ derived from 
the model of Ref. 1. Apart from the small differences, 
such as values of acoustic velocities, likely caused by our 
modification of the force constants, there is one qualita- 
tive disparity. The dispersion of the transverse acoustic 
modes in Figure 3 is clearly seen to be parabolic in con- 
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FIG. 3: The low- frequency part of calculated nanotube vibra- 
tional spectra: (a) (10,10) tube, (b) (10,0) tube. Note that 
here wave vectors fcy are measured with respect to closest V 
points of hexagons of the reciprocal graphene lattice. 



trast to a linear dispersion discussed in Refs. 1,2,5. The 
parabolic dispersion is seen over a wider range of on 
going to tubes of smaller radii, compare Figures 3 (a) 
and (b). In Sec. V we give an analytic confirmation of 
this observation within a framework of a continuum me- 
chanics. 



V. CONTINUUM ANALYSIS 

Analysis of deformations and vibrations of thin-walled 
elastic cylinders goes back as far as to Rayleigh and Love; 
see, e.g., Refs. 18-21, references therein, and Refs. 11,22 
for a dedicated analysis of vibrations. 

In the case of a planar, graphene, structure, the con- 
tinuum potential elastic energy can be written as 



I J dxdy [Ci{uoi: + Vy f 



(6) 



Uo = 



^ j dxdy[Di{wxx+WyyY 

D2{wly - W.j,,xWyy)\ . (7) 



Here x is a coordinate that would later become along the 
cylinder axis and y coordinate along the cylinder circum- 
ference, p is the mass density. Displacements fields u, v 
and w would become, respectively, parallel to the cylin- 
der axis, parallel to its circumference, and perpendicular 
to the cylindrical surface; x and y subindices denote the 
differentiation over corresponding coordinates. This type 
of deformation energy is well known in the continuum me- 
chanics of plates^^'^" and can be readily derived from the 
discrete form (1,2). In the former picture, elastic con- 
stants Ci, C2, Di and D2 in Eqs. (6,7) are expressed in 
terms of stretching and bending rigidities and Poisson's 
ratio. In the latter derivation, they would be expressed 
through constants K's in Eqs. (1,2). Note, however, that 
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FIG. 4: Branches of the vibrational spectra obtained in the 
continuum model described in the text at (a) ky = 0, (b) 
ky = l/R, (c) ky = n/R with n > 2. Calculations were 
performed for R — 6.78 A that one would have for the (10,10) 
tube. Elastic parameters used in calculations for this Figure 
would correspond to graphene velocities cia ~ 24 km/s, Cta = 
14 km/s. Dispersion of the out-of-plane graphene vibrations 
u) = dk^ was taken with S — 6 x 10~^ m^/s (see Ref. 4), and 
D2 /Di — 4C2 / Ci . Graphene results are shown by the dashed 
lines for comparison. Dispersion of the graphene transverse 
acoustic mode practically coincides with that of the nanotube 
twisting mode. 



Eqs. (1,2) should then be transformed only in the context 
of purely acoustic deformations, in which two carbons 
in a hexagon move "in-phase" . That is why we will be 
dealing here with the 3x3 matrix in Eq. (9) rather than 
with a 6 X 6 dynamical matrix. One can develop a con- 
tinuum model that would include optical "out-of-phase" 
deformations as well. The deformation energy of type of 
Eqs. (6,7) was already proved to be useful in studies of 
large deformations of carbon nanotubes.^^ 

The modification of Eqs. (6,7) upon formation of a 
cylindrical body in continuum mechanics is, e.g., dis- 
cussed in Refs. 11,20,21. As is also shown in Appendix 
A for a discrete model, this corresponds to a simple sub- 



stitution in (6,7): Vy 



w/R and Wy 



Wy 



v/R, 



R being the cylinder radius. Using this substitution, one 
can easily study the problem of small vibrations of a con- 
tinuum cylinder. The vibrational frequencies u) for the 
plane waves with a two-dimensional wave vector {kx,ky) 
are determined from the eigenvalue equation 



Md, 



(8) 



where displacement vector d = {u,v,w) and matrix 
M{kx, ky) is 



'Cikl + C2kl {Ci-C2)kxky 



{Ci-C2)kxky 
~iAkx/R 



Cikl 



C2ki 
iBky/R 



iAk^/R 
-iBky/R 
^D,{kl + klY 



(9) 

-Di/R^, C2 = C2 + D2/R^, A = 2C2-C1, 



Here Ci = Ci - 
B 

In the case of the planar structure, R 00, Eqs. (8,9) 
lead to two acoustic waves of in-plane vibrations (with 



Ci + {Di + D2/2)kl + Dikl and ^ Ci/R^. 



longitudinal cia = and transverse Cta = C^"^ veloc- 
ities) and an acoustic wave of the out-of-plane vibrations 

with a parabolic spectrum a; ~ Z)y^(fc^ -I- A:^), see Figure 
4(a). 

For a cylinder/tube of a finite radius i?, the trans- 
verse quantization imposes a restriction on values of 
ky — n/R^ where n is an integer. Of special interest 
to us here are values of ky equal to and to ±l/i?. 
The case of ky — (Fig. 4 (a)) yields two of the four 
acoustic modes of the tube with vanishing frequencies: 
(i) the longitudinal mode with a low-frequency disper- 
sion coinciding with the longitudinal acoustic wave of 
the graphene and (ii) the twisting mode whose disper- 
sion is somewhat modified from the transverse acous- 
tic mode of graphene by virtue of the elastic constant 

-1/2 

-D2: Cta = C2 ■ The longitudinal mode in the case 
of the cylinder couples with the breathing mode, whose 
frequency at fc^: = is cijf, and whose dispersion is deter- 
mined by constant Di. As a result of this coupling, an 
anti-crossing behavior of the branches arises as is clearly 
seen in the Figure, (iii) The two other, degenerate, 
acoustic modes of the tube spectrum, usually referred 
to as transverse acoustic modes for carbon nanotubes^ 
or as flexural modes in other applications, ^'^^■^'''^^ cor- 
respond to ky = ±l/i? (Fig. 4 (b)). It is apparent 
from the Figure that these modes have a parabolic spec- 
trum. In fact, one can easily show this analytically by 
the perturbation analysis of (8,9) in kx'- contributions to 
the linear coefficient in dispersion uj{kx) exactly cancel. 
The low-frequency parabolic dispersion of these modes, 
although modified by the curvature rigidity, Eq. (7), 
is mainly determined by the in-plane stretching rigid- 
ity, Eq. (6). Neglecting Eq. (7), one would obtain the 
long wavelength, for kxR ^ 1, dispersion of this mode 
as u{kx) = {2{Ci - C2)C2lCif'^ klR. Evidently, the 
parabolic character of the dispersion becomes even more 
apparent for smaller-radius tubes. These modes corre- 
spond to bending vibrations of the tube as a whole and, 
in this sense, are similar to the bending vibrations of 
rods, whose generic long- wavelength parabolic dispersion 
is well known. Similarities between Figures 4 and 3 (a) 
are obvious. Note, however, some quantitative differences 
caused by different values of effective parameters. 



VI. SUMMARY AND DISCUSSION 

We have described a simple empirical model, in which 
vibrations of the graphene and individual single wall car- 
bon nanotubes are treated and calculated on the same 
footing. Differently from the force constant model, our 
model uses only "invariant" quantities: variations of 
bond lengths, interbond and dihedral angles. As a result, 
the isometric mapping from the planar graphene onto a 
cylindrical surface of nanotubes automatically preserves 
all the right relationships between equivalent force con- 
stants. Importantly, all calculated vibrational spectra of 
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nanotubes correctly exhibit four types of acoustic exci- 
tations with vanishing frequencies (one longitudinal, two 
transverse and one twisting) , which are derived naturally 
and follow from the symmetries of the underlying system. 
Although the results obtained are largely similar to the 
earlier published, we have also found an important qual- 
itative difference. The long-wavelength dispersion of the 
transverse acoustic modes is shown to be parabolic rather 
than linear. One consequence of this is that the vibra- 
tional density of states should exhibit a one-dimensional 
singularity near zero frequency. We cannot exclude that 
there can be other physical implications of our finding, 
although this apparently is not the case for the quan- 
tized ballistic thermal conductance,^'^ for which the ex- 
act dispersion law is irrelevant. These long- wavelength 
(wavelength much larger than the tube radius) transverse 
modes correspond to bending oscillations of a nanotube 
as a whole and, therefore, are similar to bending vibra- 
tions of rods, whose low-frequency spectrum is known 
to be parabolic. The parabolic dispersion is, e.g., evi- 
dent in calculations of the lowest flexural modes of quan- 
tum wires. The origin of the parabolic dependence for 
nanotubes has been analytically illustrated using a con- 
tinuum elastic model similar to used in the analysis of 
vibrations of cylindrical shells. The continuum model 
also clarified the coupling between longitudinal acoustic 
and breathing modes resulting in the anti-crossing behav- 
ior clearly seen in the calculated spectra of nanotubes. 

The simplicity of our model allows us to easily calcu- 
late vibrational spectra of nanotubes of arbitrary chiral- 
ities. To better handle such spectra, we employed an 
alternative definition of the nanotube unit cell with only 
two carbons per cell. This definition reveals a "true" 
longitudinal periodicity of carbon nanotubes that can 
be substantially shorter than used in the conventional 
definition.^ 

Being a result of the straightforward isometric map- 
ping, the model presented in this paper, has its quali- 
tative limitations. For instance, it does not automati- 
cally yield the curvature- induced softening of the phonon 
modes observed in ab initio studies such as that of 
the twisting mode,^ or more complex dependencies for 
the breathing mode.^'^'* Of course, these effects are not 
deducible from the graphene properties. To describe 
them, the model would need to include both the relax- 
ation of nanotube geometries (mapping would not be 
exactly isometric) as well as an explicit dependence of 
the elastic energy on nanotube radius and with the axial 
anisotropy.^^ We discussed such a generalization in the 
context of uniform deformations,^^ which, e.g., leads to 
the chirality-dependcnt stiffness of nanotubes.^® We be- 
lieve the present model can also be further developed in 
this regard and correspondingly parameterized by com- 
parison with ab initio calculations. A more accurate pa- 
rameterization is, needless to say, required even for the 
present model. Such quantitative aspects have not been 
pursued in this paper. 
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APPENDIX A: A RING EXAMPLE 

Here we study in more detail a simpler illustrative ex- 
ample of the relationship of vibrations of an infinite linear 
atomic chain and a ring of atoms. For clarity, displace- 
ments of atoms are restricted to the plane in which a ring 
and chain belong. Atom-atom interactions correspond 
to formation of bonds and result both in stretching and 
bending rigidity of a linear chain. Figure 5 is a picture 
of the ring consisting of TV = 10 atoms, which is assumed 
to preserve the nature of interactions in the chain. The 
figure shows a common system of coordinates (y, z) as 
well as local systems {v, w) related to the tangential and 
normal displacements of the corresponding atom. Obvi- 
ously, the kinetic energy 

n n 

where n is the 1-d positional index of the atom along the 
ring circumference. We first write the potential energy 
of a linear chain in an invariant form as 

U^Y^Sll + ^Y^^xl (A2) 

n n 

where the first term describes the stretching and second 
the bending rigidity. Then we assume that the same 
functional form holds for the ring. From the geometry of 




FIG. 5: A ring of A^ = 10 atoms connected by interatomic 
"bonds" (thicker lines). The stretching rigidity corresponds 
to variations of the bond lengths I, the bending rigidity to 
variations of the inter-bond angles x- Number of atoms N 
(or the radius of the ring) determines the angle (fi = 2n/N. 
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FIG. 6: The vibrational spectra of the linear chain (solid lines) 
and of the ring (filled circles) of atoms. The parameters used 
for the plot are K/M = 1 and K^/M = 0.1. 



and the results are shown in Figure 6 with filled 
circles. In accordance with Eq. (A6), it is now a set 
of discrete frequencies. The coupling between normal 
and tangential displacements resulted in important 
qualitative modifications of the spectrum, which have 
been obtained exactly and naturally. Particularly, one 
notices 3 zero-frequency modes. Two of them (with i = 1 
and i = A'^— 1, ori = — 1, in Eq. (A6) - precisely one 
wavelength on the ring circumference) correspond to the 
displacements of the ring as a whole in two orthogonal 
directions. The other (with i = 0) is a pure tangential 
mode describing the rotation of the ring as a whole. The 
k = Q finite-frequency mode, on the other hand, is a pure 
normal, breathing, mode that "borrowed" its strength 
from the parent stretching oscillations of the linear chain. 



Figure 5, the variation of the bond lengths, as a function 
of the local displacements, is 

Sin = (wn+i-i'«)cos((/3/2) + (w„+i+w„)sin((^/2), (A3) 

while the variation of the intcr-bond angles 

I ■ SXn = (2W„ - Wn+l - Wn-l) C0s((^/2) 

+ K+i - Wn-i)sm((^/2). (A4) 

For the linear chain {(p = and w„ = ?/„, Wn = Zn), 
Eq. (A3) describes a conventional longitudinal stretching 
while Eq. (A4) would yield a conventional local curva- 
ture. Tangential and normal vibrations in a linear chain 
arc fully decoupled. The corresponding two branches of 
the vibrational spectrum are simply 

c.2(fc) = (4K/M)sin2(fca/2), 

Lolik) = (mKb/M) sm\ka/2), ^ ^ 

where a = / is the distance between the atoms along the 
chain. This spectrum is shown in Figure 6 with solid 
lines. 

In the curved system (our ring with ip ^ 0), on the 
other hand, tangential and normal displacements are 
coupled, as is well known in the elasticity theory for 
curved surfaces"'^^'^^ and clearly seen in Eqs. (A3,A4). If 
we were to use the common system of coordinates for dis- 
placements {um Zn), the atom contributions to the poten- 
tial energy, Eq. (A2), would be position dependent and 
the translational invariance with respect to n n + 1 
would be lost. With the local coordinates {vn,Wn), this 
invariance is preserved and one can directly use the con- 
ventional transition to A:— states, which would now be 
quantized as 

k^{2T:/Na)i, i = 0, ...iV-1. (A6) 

The unit cell length a > I here is now the distance 
between the atoms of the ring along its circumference. 
The derivation of the spectrum from the equations of 
motion follows straightforwardly from Eqs. (A1-A4), 



Equations (A3,A4) can be used to study the continuum 
limit, when, keeping the same radius R, we increase the 
number of atoms N oo and bond length 1—^0. Then, 
evidently, 

Sl/l^[dv/dy + w/R], 

Sx/l ^ —dldw/dy — v/R] /dy. ^ ' 

The expressions in brackets in Eq. (A7) provide a recipe 
for a transition from the continuum model of a chain to 
the continuum model of a ring. The subsequent analysis 
of vibrations is straightforward with the wavevector k 
quantized as kR = n (n being an integer) yielding three 
zero-frequency modes and the breathing mode as in the 
discrete case above. 

It should be noted here that in reality the very values 
of the unit cell length a and elastic constants K , Kb can 
in fact somewhat differ for the ring and the linear chain, 
and the difference would be iV— dependent as determined 
by the equilibrium bond length in the ring. In this sense, 
what is compared in Figure 6 is the structure of the spec- 
tra for the same values of the essential parameters upon 
isometric mapping of the linear system onto a circular 
ring. In addition, the elastic energy of the ring could in 
general contain terms absent in the energy of the chain, 
such as SxnSln- They are not deducible from the func- 
tional form of the linear system and would have to be 
explored on their own. We also note that Eq. (A2) can 
be generalized to include longer range interactions. 

APPENDIX B: TRANSVERSE QUANTIZATION 
AND UNIT CELLS 

Two carbons connected by Ch (Eq. (3)) on the 
graphene plane correspond to the same carbon on the 
nanotube after wrapping. The effective 2-d cyclic condi- 
tion can therefore be written as 

R-nm + C/i = R„m, R„m = nai -I- ma2, (Bl) 

R-nm being the position of one of the carbons (there are 
two of them) of an arbitrary unit cell of graphene. If k 
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FIG. 7: Unwrapped (3,2) tube with carbons rearranged to 
form a horizontal strip here. The tube axis is shown by the 
dash-dotted line. Here dc = 1 and the T vector, Eq. (B2), 
is defined by the {P,Q) pair (2,f). Small circles connected 
by thinner lines indicate a sequence of carbons reached by 
the T— translation. It is transparent that all hexagons of the 
strip would be visited this way, forming an effectively 1-d 
enumeration of Eq. (B5). The conventional vector Td ~ 
7ai — 8a2 and such a unit cell would have 38 hexagons. 

is the 2-d wave vector of the band excitations, then in 
the infinite graphcne plane it would have two indepen- 
dent continuous components. In nanotubes, the cychc 
condition leads to the transverse quantization of Eq. (5). 
In other words, wave vectors k = -I- ky allowed by 
Eq. (5) lie only on certain quantization lines in the recip- 
rocal plane of graphene, kj^ being quantized according 
to Eq. (5) and k|| being actually a 1-d continuous wave 
vector parallel to the tube axis. 
Let us define vector 

T = Pai + Qa2 (B2) 

with integer P and Q. From Eqs. (3) and (B2), one finds 

ai = i-QCh + A/T)/A, el2 = (PC,, - iVT)/A, (B3) 

where 

A =: MP - NQ. (B4) 

Evidently, if one can find such P and Q (of course, we are 
interested in the "smallest" P and Q) for a given tube 
{N,M) that A = ±1 in Eq. (B4), then original primi- 
tive vectors ai and a2 in Eq. (B3) will be represented 
through integer amounts of C,, and T. Correspondingly, 
an arbitrary vector R„m, Eq. (Bl), will be expressed 
through integer quantities of Ch and T as well. Since the 
cyclic condition defines Rnm in Eq. (Bl) with accuracy 
to Ch, this would actually mean that the position of any 
hexagon in the unwrapped tube is determined through a 
single vector T: 

R„„, ^ R, = iT, (B5) 

where, e.g., i — nM — mN for A = 1. One can think of 
the corresponding unit cell built of Ch and T that would 
contain only 2 carbons in the cell. Of course, vector 
T does not have to he parallel to the tube axis. The 




FIG. 8: Unwrapped (4,2) tube with carbons rearranged to 
form a horizontal strip here. The tube axis is shown by the 
dash-dotted line. Here dc = 2 and the T vector, Eq. (B2), is 
defined by the "reduced" {P,Q) pair (1,0). Two sequences of 
carbons resulting from T— translations are shown by thinner 
lines: small circles of one sequence are connected by the solid 
line, and circles of the other by the dashed line. Translation 
from one sequence to the other is achieved with vector Ch/2. 
All hexagons of the strip would be visited this way, forming 
effective enumeration of Eq. (B8). The conventional vector 
Td = 4ai — 5a2 and such a unit cell would have 28 hexagons. 

projection of T on the nanotube axis Til = V3a'^A/2Ch 
- it would be directed in opposite ways for A = 1 vs 
A = — 1 - and its modulus determines the corresponding 
longitudinal period 

Oil = ^ia^l^Ch, (B6) 

where Ch = \Ch\- Note that vector T is different from 
the symmetry vector defined in Ref. 1. 

It is easy to see that the picture described in the pre- 
vious paragraph is indeed realized whenever dc = 1 - we 
will call it the irreducible case - where dc is the greatest 
common divisor of TV and M . (Considering only posi- 
tive N and M does not restrict the generality.) Then, 
in fact, any integer value of A in Eq. (B4) can be estab- 
lished with an appropriate choice of P and Q. Examples 
of positive (P, Q) pairs satisfying A = 1 are listed here 
as (7V,M) ^ (P,Q): (3,1) -> (1,0), (3,2) -> (2,1), 
(5,2) (3,1), and (10,9) (9,8). An illustration for 
the (3,2) tube is shown in Figure 7. A visual picture of 
the irreducible case is that of a 1-d chain (of period |T|) 
that is wrapped around the nanotube cylinder with an 
appropriate helix angle, as is clearly seen from Figure 7. 

"Reducible" cases - with dc > 1 - can be described in a 
similar but somewhat more involved fashion. One would 
extract the irreducible structure factors A^i and Mi: 

N^dcNi, M^dcMi, 

so that the greatest common divisor of A^i and Mi equals 
1 . Then the procedure described above for the irreducible 
case can be applied for the geometry {Ni,Mi) resulting 
in the first translational vector T. In other words, inte- 
gers P and Q of Eq. (B2) should satisfy condition 

MP - NQ = dc. (B7) 

One however could not visit all hexagons by using only 
so defined T. The needed second translational vector 
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FIG. 9: An illustration of the range of variation of ky for 
(3,2) and (4,1) tubes. In the former case K2 = 2bi — 3b2, 
in the latter case K2 = bi — 4b2. The corresponding ranges 
are shown as thick lines between the centers of hexagons of 
the reciprocal lattice of graphene. Each line crosses several 
hexagons. If all inside-a-hexagon segments are displaced to 
one hexagon with a graphene reciprocal vector, they would 
form a traditional picture of quantization lines in the first BZ 
of graphene. 

can be found as Ch/dc- The resulting enumeration of all 
hexagons will then read as 

R,„„ ^ Ry = iT + jCh/d,, J = 0, 1, ... 4 - 1. (B8) 

Index i here is, as in Eq. (B5), an arbitrary integer 
that would define unique carbons while j results only 
in dc unique translations due to the cyclic condition in 
Eq. (Bl). An illustration for the (4,2) tube is shown in 
Figure 8. A visual picture of the irreducible case is then 
of dc 1-d chains (of period |T|) that are wrapped around 
the nanotube cylinder with an appropriate helix angle. 
The elementary unit cell contains two carbons and built 
of vectors T and Ch/dc- Longitudinal periodicity will 
again be determined by the projection of vector T on 
the nanotube axis which is Eq. (4) becoming Eq. (B6) at 
4 = 1. 

For the reciprocal lattice vectors, defined through Ki • 
C,, = 27r, K2 • T = 27r, Ki • T = 0, and K2 • C,, = 0, one 
readily obtains 

Ki -Qbi+Pb2, K2 = Mibi - iViba, 

where bi and b2 are primitive graphene reciprocal vec- 
tors defined with respect to ai and a2. One easily finds 
that IK2I indeed coincides with 27r/a||. An illustration 
for cases of (3,2) and (4,1) tubes is given in Figure 9. 

This way one arrives at the band excitation picture 
with 2dc (for each component of the wave function or 
polarization) continuous bands corresponding to 2 atoms 
in the unit cell and dc transverse quantization levels, e.g. 
with quantum numbers 

/ = -4/2 + l,...,0,...4/2. (B9) 



The latter define the quantized component of the parent 
2-d quasi-momentum k perpendicular to the tube axis: 

fci = i2n/Ch)l. (BIO) 

The continuous component fc|| parallel to the tube axis, 
on the other hand, would be defined within a BZ whose 
width is 

2Tr / a\\ inCh/VSdca^. (Bll) 
One can easily calculate that the resulted number of 
continuous bands and the width of the BZ lead to the 
correct total number of states in the system, which 
corresponds to 2Ch/V^a^ hexagons per unit length of 
the nanotube. 

The conclusion that there are only dc "truly unique" 
quantization levels can be confirmed another way as well. 
If G = nibi + n2b2 is the vector of the reciprocal lattice 
of graphene, then ki = k + G is physically equivalent 
to k. In general, k and ki may correspond to different 
quantum numbers in Eq. (5), say / and h. One derives 
the integer difference of these quantum numbers as 

dl^l-li^ mN + n2M = dc{niNi + 7i2Afi). (B12) 

Once again, the last, "irreducible", factor in Eq. (B12) 
can take any integer values: 

mNi + nzMi = 0, ±1, ±2, . . . , 

with appropriate choices for integer ni and 712 . It is then 
clear that the number of unique sub-bands in the ex- 
tended scheme equals precisely dc ~ one can, e.g., choose 
quantum numbers of Eq. (B9) for "independent" sub- 
band indexing, all other values of I would be reducible 
to independent values with the appropriate adjustment 
for the one-dimensional continuous quasi-momentum k\\ . 

It is worth mentioning that in the {N, M) family of 
nanotubes with iV > Af , it is the armchair and zigzag 
nanotubes that have maximal number of transverse 
quantization levels dc — N. The tubes immediately 
next to them, {N,N — 1) and {N, 1), on the other hand, 
would have only one quantization level. The armchair 
tubes have then the longitudinal period of a/2, twice as 
small as is with the conventional definition.^ One can 
easily see it with T vector defined by the (P, Q) pair (1,0). 
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